Magnetic and transport properties of the one-dimensional ferromagnetic Kondo 

lattice model with an impurity 
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We have studied the ferromagnetic Kondo lattice model (FKLM) with an Anderson impurity on 
finite chains with numerical techniques. We are particularly interested in the metallic ferromagnetic 
phase of the FKLM. This model could describe either a quantum dot coupled to one-dimensional 
ferromagnetic leads made with manganites or a substitutional transition metal impurity in a MnO 
chain. We determined the region in parameter space where the impurity is empty, half-filled or 
doubly-occupied and hence where it is magnetic or nonmagnetic. The most important result is 
that we found, for a wide range of impurity parameters and electron densities where the impurity 
is magnetic, a singlet phase located between two saturated ferromagnetic phases which correspond 
approximately to the empty and double-occupied impurity states. Transport properties behave 
in general as expected as a function of the impurity occupancy and they provide a test for a 
recently developed numerical approach to compute the conductance. The results obtained could be 
in principle reproduced experimentally in already existent related nanoscopic devices or in impurity 
doped MnO nanotubes. 

PACS numbers: 75.30.Mb,75.47.Lx,75.47.-m,75.40.Mg 



I. INTRODUCTION 

Manganese oxides, such as Lai-^Ca^MnOs, commonly 
referred to as manganites, have attracted an intensive 
theoretical and experimental effort, 1,2 mainly due to 
their property of colossal magnetoresistance^ and its 
consequent applications to magnetic recording devices. 
General applications of the ferromagnetic (FM) metallic 
phase of manganites belong to the field of spintronics^ 
where the spin of the electrons is exploited in addition to 
its charge. 

A simple spintronics device which is relevant for the 
present study is formed by a quantum dot (QD)j* a 
nanometer-scale box, connected to two FM leads' This 
device can act as a spin valve^ or a spin filter. Ferro- 
magnetic metals (Co, Pd-Ni) or diluted magnetic semi- 
conductors, such as GaMnAs, are employed as leads. 
Alternatively, manganites are also used as FM leads in 
spintronic o °i because of its high polarization. 

Manganites are usually described by the ferromagnetic 
Kondo lattice model (FKLM) in which the conduction 
sites represent the orbitals e g and the localized spins 
the orbitals t2 g ~ The QD will be described as a sin- 
gle Anderson impurity. The spin valve with manganites 
as leads corresponds then to a FKLM with an Ander- 
son impurity which is the model we will study in the 
present work. Moreover, we will consider this model in 
a one-dimensional (ID) space. One should keep in mind 
also that the magnetoresistance of manganites is usu- 
ally applied in multilayer heterostructures FM/M/FM or 
FM/I/FM (M: metal, I: insulator) which can be consid- 
ered as ID systems in the direction perpendicular to the 
interfacefi 2 - 

This model, in addition to its application to a wide 
class of devices, could describe other more conventional 



condensed matter systems such as a transition metal ion 
such as Cu replacing Mn in a manganese oxide chain.— 
It is well-known that a single impurity could lead to in- 
teresting and important local or short-range effects in 
magnetic systems in low dimensions^ These effects can 
in turn modify the long-range physics of such systems for 
a finite density of impurities. 

The main purpose of this work is to search in the pa- 
rameter space of the model for phases in which the satu- 
rated FM is reduced to a partially polarized FM, or even 
to a nonmagnetic state, upon the introduction of an An- 
derson impurity. This problem would be the analog for a 
ferromagnetic chain of the effect that causes a magnetic 
impurity in a paramagnetic metallic chain, that is, the 
paradigmatic Kondo effect^ We would like to empha- 
size that finite size effects due to our finite-cluster calcu- 
lations could be relevant both to describe mesoscopic de- 
vices and to capture local or short-range features caused 
by an impurity in a Mn-0 chain. For completeness, since 
the model studied may be applied to electronic devices, 
we will compute the conductance through the QD but 
clearly the study of transport properties is not the main 
motivation for the present work. In any case, even though 
the physics found for most of the parameter space corre- 
sponds to the saturated FM phase and hence transport 
properties can be recovered by a spinless fermion model, 
the impurity-FKLM is an interesting testing ground for 
the quite recent numerical techniques we will employ for 
this study. 



II. MODEL AND METHODS 

Hence, in this article we will study a one-dimensional 
FKLM with an Anderson impurity located in the center 
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FIG. 1: (Color online) Picture of model (JTJ. 

of the chain (see Fig. [1} . Then, the model is defined by 
the Hamiltonian: 

H = - t ^2 {c\ +la Cia + H.c.) - J H ^ S z • S; 

»>Q,i<— l,ir Z#0 

+ C/7i ,Tn 0: x + en (1) 

where the notation is standard. The Anderson impurity 
or "QD" , with parameters U, e, is located at site "0" and 
is connected to the rest of the system with a hopping 
t'. The "leads" (i ^ 0) correspond to the FKLM with 
the Hund's rule exchange coupling Jh > 0. Si is the 
spin operator for the localized spin-1/2 orbital and Si 
is the one for the conduction electron at site I (I ^ 0). 
to = 1 is adopted as unit of energy, and we take t' = 0.4 
throughout. Model Q} will be termed "QD-FKL" model. 

The pure single-orbital FKLM or Kubo-Ohata model^ 
has been extensively studied, particularly using numeri- 
cal techniques^ and its phase diagram for various spa- 
tial dimensions and values S of the localized spins has 
been determined.— Even in the simplest case of ID and 
spin-1/2 localized spins, the model reproduces qualita- 
tively the main features of manganites. In the following 
we will work in the metallic FM phase of the FKLM, 
typically, the density of conduction electrons n < 0.6, 
Jh = 20 (all coupling constants are expressed in units 
of to). The on-site potential e and Coulomb repulsion U 
are the main variables whose effects we want to study. 
In a heterostructure e would be fixed by chemistry but 
in a spin valve it would correspond to the gate voltage 
which can be varied at will. In order to detect any de- 
parture from the fully polarized FM state it is essential 
to work in the subspace of total S z = (1/2) for even 
(odd) number of electrons. 

We denote with L the total length of the system in- 
cluding the impurity site. Open boundary conditions 
(OBC) were adopted in the lattices studied except oth- 
erwise stated. Small clusters with L up to 12 will be 
studied using exact diagonalization (ED) with the Lanc- 
zos algorithm. Larger clusters will be solved using den- 
sity matrix-renormalization group (DMRG)i^ For calcu- 
lations in the subspace of maximum total z-component 
of the spin, 5 Z , i.e., saturated ferromagnetism, we used 
completely independent ED and DMRG codes for the 
spinless fermion model. The ED and DMRG codes for 
the FKLM were thoroughly checked, in the first place by 
reproducing results in the literature. In the second place, 
by comparing results obtained by both techniques in 
small clusters and finally, by comparing results for large 
chains between FKLM and the spinless fermion model in 



the case of maximum S z . Here we would like to stress the 
fact that convergence to the ground state, both with ED 
and in the diagonalization of the superblock Hamiltonian 
at each iteration of DMRG is extremely slow. This is al- 
ready known for ED studies of the Hubbard model with 
very large U jt and small doping, that is in the proximity 
of the Nagaoka phase. DMRG studies for the Kondo lat- 
tice model (KLM), with both ferro- and antiferro- mag- 
netic exchange coupling, have been in general restricted 
to smaller chains than for the Hubbard model. In fact, 
even in ID, the DMRG treatment of KLM has the level 
of difficulty of an interacting system on a two- leg ladder. 
The convergence is even worst for the FKLM where previ- 
ous studies have been limited to L k 36 with a discarded 
weight of order 10 _5 A 9 - Last but not least, the presence 
of impurities makes the convergence more difficult par- 
ticularly for DMRG. In our calculations, with a retained 
number of M » 400, the truncation error is negligible 
(w 10~ 14 ) for L w 20 in the regions close to saturated 
ferromagnetism but it drops to 10~ 10 in the nonmag- 
netic region. In the case of the spinless fermion model, 
for L w 20 and M w 400, the precision in energy is at 
least 12 digits. In the parameter regions where total spin 
S takes its maximum possible value 5 mal , the energy 
in the S z = subspace reproduces the value obtained 
in the S z = S ma x subspace using the spinless fermion 
model within at least 9 digits. In any case, the limited 
precision within DMRG depends essentially on the lack 
of convergence in the diagonalization of the Hamiltonian. 

The conductance will be estimated by a numerical 
setup 20 ' 21 in which a small bias voltage is applied to 
the left (L) and right (R) leads, with AV = Vr - Vl, 
(Vr = —Vl), at time t = 0. The current J{t) induced 
by this voltage on each bond connecting the QD to the 
leads is computed with the time evolution formalism both 
within ED or DMRGi^ 2 . This numerical setup is equiva- 
lent to the systems which were treated analytically using 
the Keldysh Green functions formalism. 2 - These results 
for out-of equilibrium, with interacting QD, contain as 
particular cases the ones for noninteracting systems de- 
scribed by the Landauer formula* 2 ^ These analytical re- 
sults, both for interacting and noninteracting QDs, were 
recovered using this numerical setup and time-dependent 
DMRG^ffi The advantage of this numerical proce- 
dure is that it can be extended with no formal limitations 
to study the case of interacting lead^ as will be done in 
the present work. 

In principle, one could adopt as a measure of the 
conductance the maximum of J(t)/AV. It has been 
shown that this recipe provides correct results for the 
conductance when the maximum of J(t) corresponds 
to a "plateau" which appears in high-precision calcu- 
lations using the "adaptive" timc-dependent-DMRG on 
large clusters. 20 In the following we adopt the less precise 
"static" algorithm 25 which still gives qualitatively correct 
results, particularly if relatively small clusters are consid- 
ered, but is much faster than the "adaptive" scheme, thus 
allowing to explore a wider range of couplings and den- 
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FIG. 2: (Color online) (a) QD occupancy as a function of e, at 
several electron densities n. (b) AE — E(S Z ) — E(0), versus 
S z , n = 0.5. (c) J(t)/AV for n = 0.4, e = -4.5, -5.5 and 
— 12 (full curves). Sinusoidal fits are shown with dashed lines. 
Curves for e = —4.5 and —12 have been multiplied by 2 and 
20 respectively, (d) Conductance, as a function of e, for the 
same electron densities as in (a). L = 10, Jh = 20, t' — 0.4, 
e = -17/2. 



sities. In the case of ED, the time evolution is exactly 
computed in the full Hilbcrt space of the system. The 
time-evolution of the ground state is given by ^(i+r)) — 
e- mr \*(t)), where H = H + V L N L + V R N R , N L , N R 
are the electron occupancies of the left and right leads 
respectively, and |*(t = 0)} = |* ), ^ol^o) = ^o|*o)- 
\ty(t + r)) was computed using the Krylov algorithm. 27 
All the results reported below correspond to AV = 0.01, 
r = 0.1. 

The dynamical impurity magnetic susceptibility and 
dynamical magnetic structure factor (defined for con- 
venience in Section IIV|) are computed within ED and 
DMRG using the standard continued fraction formalism. 
In the case of DMRG we again choose the "static" formu- 
lation which although less precise is enough to determine 
the presence or absence of a peak at the bottom of the 
spectrum. 

We would like to stress that the most important results 
reported in this article correspond to static properties, 
that is, ground state energies and spin-spin correlations, 
where the precision of DMRG is maximal. For all quan- 
tities studied, the results obtained with ED are precise 
to precision machine. 



III. RESULTS AT THE SYMMETRIC POINT 

Let us start to analyze results for the QD-FKL model 
in the L = 10 cluster obtained by ED. We consider 
in the first place the case of the "symmetric point" , 
e = —U/2. The symmetric point of an Anderson im- 
purity is in principle the obvious place to look for a mag- 
netic impurity placed in a noninteracting chain. How- 
ever, this is not the case for the present model. It 



can be seen in Fig. [UJa) that the QD or impurity oc- 
cupancy, and hence its actual magnetic or nonmagnetic 
character, experiences a sharp crossover as a function 
of the on-site potential of the impurity. For values of 
e larger than a crossover value e*, tiqd « 0, and for 
values of e < e* , tiqd « 1. e* may be defined as the 
value of e at which uqd = 0.5. This crossover can 
be understood by examining two variational states in 
the atomic limit. One, with energy E\ = —JuTie/A, 
(n e is the number of conduction electrons) where all 
electrons are located on the leads and ferromagnetically 
aligned with the localized spins, and the other with en- 
ergy E 2 = — Jff(n e — l)/4+e, where one electron has been 
moved from the leads to the QD. The crossover between 
both variational states at e* ar — —Jh/^ is quite close to 
e* as shown in Fig. HJa). The dependence of e* with n 
is mainly due to the kinetic energy which can be easily 
computed within the spinlcss fcrmion model to which the 
FKLM is reduced in the saturated FM state, i.e. when 
total spin S = S max — S* lax — n e /2. In this case, ne- 
glecting the term with t', £* spin i ess is equal to the single 
particle energy of the top of the band which increases 
with n and is exactly zero at n = 0.5. The connection 
between both models implies that e S piniess — £ + Jh /4. 

It is important to notice at this point that although 
the pure system is in the saturated FM state for the 
densities studied, for some impurity parameters the im- 
purity may drive the system into partially polarized FM 
states with total spin S < S max . In fact, as shown in 
Fig. Hfb) for n — 0.5, there is a region of e, close to e*, 
where S = (1/2) for even (odd) n e . We have observed 
this nonmagnetic state for other chain lengths and den- 
sities. Although the difference in energy between states 
with different S z is very small, these results strongly sug- 
gest that there are regions in parameter space where the 
impurity causes a breakdown of the fully saturated FM 
state. This possibility will be thoroughly examined in 
the next Section. 

Let us discuss in detail how the conductance G is de- 
termined. In Fig.[2{c), it is shown J(t)/AV (J(t) is the 
average of the current on the two bonds connecting the 
QD to the leads) which presents the typical oscillatory 
behavior. This oscillatory behavior follows from the ex- 
pansion of e _lWr in eigenvectors of H, which for small 
Ay are adiabatically related to those of Hq. We would 
like to emphasize that results depicted in Fig. \T[.c) are 
exact, i.e., no truncation of the Hilbert space was per- 
formed. Then, we fit each curve by a sinusoidal and we 
adopt G as the amplitude of this sinusoidal. In this small 
cluster, but also for L = 20, a single sinusoidal gives a 
reasonable fitting of J(t)/AV for most of the cases stud- 
ied, particularly near e*. Although this procedure is not 
very precise, it gives qualitatively correct results as we 
discuss in the following. 

Results for the conductance are shown in Fig. [U(d) for 
various densities as a function of e. G is only different 
from zero at the crossover between the region of empty 
QD (e >> e*) and the region of half-filled QD (e << e*) 
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FIG. 3: (Color online) (a) QD occupancy and (c) conduc- 
tance, as a function of e, for several electron densities n indi- 
cated on the plot, (b) J(t)/AV for n = 0.4, e = —5, —5.5 and 
—7 (full curves). Fits to a sine function are shown with dashed 
lines. Results for L = 20, J H = 20, t' = 0.4, e = -(7/2. (d) 
Spin-spin correlations along the conduction chain, for the pure 
FKLM (stars) and in the presence of an impurity for various 
values of e indicated on the plot, L = 19, n = 0.421. The 
reference site is located at the center of the chain and the 
normalization (So So) = 1 was adopted. 



and it has a sharp peak at e* with a width approximately 
equal to the bandwidth of a tight-binding model on the 
leads, 4tQ. This behavior is what one would expect for 
the spinless fermion model. The determination of the 
variation of the maximum conductance with density n, 
which would require calculations on a finer mesh in e, is 
out of the scope of the present study. 

Let us now discuss results for L = 20, obtained with 
DMRG, also at the symmetric point. In Fig. [3Ja), it can 
be seen that, as for L = 10, the impurity occupancy uqd 
experiences a sudden change as a function of the on-site 
potential e. This crossover is located approximately at 
e* = —Jji/i as argued before, even inside the incommen- 
surate (IC) phase (but not strictly at n = 1), which is 
expected since the variational states are independent of 
the underlying FM or IC order. As in the L = 10 cluster, 
the location of this crossover shifts to larger values of e 
as the density is increased. At the IC-FM crossover for 
Jh = 20, n 0.55, e* experiences a somewhat larger 
increase and then it remains relatively unchanged up to 
half-filling. In the IC region of course the kinetic energy 
is no longer approximated by the spinless fermion model. 
In fact, it is easy to realize that the kinetic energy ver- 
sus e follows an opposite behavior as the IC-FM border 
is crossed. 

The computation of the conductance follows the steps 
previously outlined. In Fig. OJb), J(t)/AV is shown for 
n = 0.4 and several values of e. It can be seen that 
in spite of the approximate nature of the computation 
of the time evolution in a truncated Hilbert space, J(t) 
is clearly well fitted by a single sinusoidal, particularly 
close to e*, and these oscillations have similar behaviors 



FIG. 4: (Color online) (a) Variational states for model Eq. (TTJ 
in the U-e plane. Vertical thick lines indicate approximately 
the region with minimum 5* for some values of U, n = 0.4, 
discussed in the text, (b) QD occupancy and (c) conductance 
as a function of e, for U — 4, and various electron fillings n. 
Results for L = 10 chain except otherwise stated. 

as a function of e* as earlier for the L = 10 chain. In 
Fig. [3tc) we show the resulting G as a function of e and 
for several densities. As for the smaller lattice L = 10, 
the conductance is different from zero only for ewe*, 
with a peak at e*<2£ 

In spite of the extremely slow convergence in DMRG 
calculations, by keeping 450 states we were able to find 
that E(S Z = 0) < E(S Z = S max ) for e = -8, thus sug- 
gesting for L = 19 a similar behavior to that shown in 
Fig. Wt c ) for L — 10. Further indications of this be- 
havior can be obtained by examining the ^-component 
of the spin-spin correlations (SjSq), where the reference 
site "0" is the center of the chain and j labels conduction 
sites. Due to the large value of Jh, the correlations be- 
tween the impurity site and the localized spins have the 
same qualitative behavior so in this and in the following 
section we will only consider the correlations between the 
impurity and the conduction sites. These correlations, 
shown in Fig.fS^d) for n = 0.421, clearly depart from the 
correlations in the pure system as |e| is increased. This 
behavior indicates that the ground state computed by 
DMRG is a mixture of states very close in energy which 
depart from the saturated FM state. Then, this behav- 
ior of (SjSq) with e, which follows the same trend as 
the one for the L = 10 chain, suggests that also for the 
L = 19 chain the ground state S < S max . Notice also 
that in this \ow-S region, (S z Sq) do not show any trace 
of antiferromagnetic order. 



IV. BREAKDOWN OF THE FERROMAGNETIC 
STATE AWAY FROM THE SYMMETRIC POINT 

Let us discuss the consequences of these results for de- 
vices where the leads are made with manganites. Since 
Jh /to in manganites has been estimated of the order or 
larger than 10^ then a value of U* — -2e* > 5 in the 
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FIG. 5: (Color online) (a) AE = E(S^ ax ) - E(S^ in ), as 
a function of e, for (7 = 4, and various electron fillings n. 
Results for L = 10 chain except otherwise stated, (b) Upper 
(solid line) and lower (dashed line) boundaries of the region 
where the ground state S — 0(1/2), for n — 0.4, as a function 
of L, for U = 4, and n = 0.4. (c) < S z > (normalized to 1) 
at the impurity site (full symbols) and total spin normalized 
to S m ax (open symbols), as a function of e, for U = 4. Results 
for L = 11, n = 0.364 (circles), L = 19, n = 0.421 (diamonds), 
and L — 20, n = 0.4 (triangles). 



QD (see Fig.^Ja)) would be required for uqd ~ 1. This 
value of U is somewhat larger than the one in materials 
employed in the QD, such as semiconductors or carbon 
nanotubes. It is necessary then that the device could be 
operated away from the symmetric point. More impor- 
tantly, a smaller U could imply a larger effective coupling 
between the impurity and the conduction sites, assuming 
that to lowest approximation the relation J e ff w t' 2 /U is 
still valid for the present model for U > t' . In support of 
this hypothesis, we have observed that the spin-spin cor- 
relation between the impurity and its nearest neighbor 
site becomes more negative with decreasing U at fixed e. 
Then, by working with a larger J e f f we could expect to 
be more able to detect the presence of the nonmagnetic 
phase which was suggested by the results found in the 
previous section. For these reasons, in the following we 
adopt a moderate value of the on-site Coulomb repul- 
sion, U — 4, and we study the properties of the model 
for variable e, i.e. following a vertical line in Fig.[3Ja). 

The electron occupancy at the QD, shown in Fig.[3Jb) 
for L = 10, presents now three regions where ngn is ap- 
proximately 0, 1, and 2 as e decreases. The crossovers 
among these regions are located near the variational es- 
timates £q i ~ — Jff/4, and e* 2 w — J#/4 — U which 
are shown in Fig. BJa). As it can be seen in Fig. BJc), 
the conductance G consistently with the results shown in 
Figs. d) and 03c) at the symmetric point, presents 
sharp peaks at the crossovers between regions with differ- 
ent tiqd, i.e. when uqd w 0.5 or 1.5. Notice that G has 
a larger value between the peaks for n = 0.4 compared 
with the one for n = 0.3, which is essentially zero. 

Figure [5] contains the most important result of our 



work. In Fig. Eta) we plot AE = E(S max ) - E(S min ), 
where S min = S mm = (1/2) for even (odd) n e , as a 
function of e in the L = 10 chain. It can be clearly seen 
that the ground state S is smaller than S max for the var- 
ious densities considered for e < €q x . In this case, AE is 
quite large and we were able to obtain for L = 20 results 
very close to those for L = 10, as shown for n = 0.4, sug- 
gesting that this feature is at least not an artifact of this 
small chain. More interesting is the fact that inside the 
tiqy) = 1 region there is an interval in e, which depends 
on the density, where S = S m i n , as shown in Fig. E^b) 
for n = 0.4 as a function of 1/L. This state appears be- 
tween two saturated FM states, 5 = S max — 1 below the 
lower boundary line because of the double occupancy of 
the QD, and S = S max above the upper boundary line 
where the QD is empty. The error bars quoted in the 
plot correspond to the grid adopted in the e axis. For 
the largest chain studied, L = 30, we have only computed 
the energies in the S z = S max , S max - 1, S max - 2 and 
Smax — 3 subspaces. Decreasing precision as L increases 
prevents us to take a finer grid close to the crossover be- 
tween different regions and hence error bars are larger. 
An extrapolation to the bulk limit would not be reliable 
with these error bars. For L = 10, n = 0.3, the region in 
e where S — 1/2 shrinks to [—9.3,-8.6]. It is tempting 
to relate this smaller interval for n = 0.3 with respect 
to the one for n = 0.4 with the behavior of G noticed 
above but further analysis would be needed to confirm 
this possibility. In Fig.^c) we show that the presence of 
a state with S — S m i n corresponds to a strong magnetic 
character of the impurity with < S z > approaching its 
maximum value of one (according to the normalization 
adopted). This feature has been observed for all lattices, 
densities, and values of U studied. The intervals in e 
where S = S m i n for various values oft/, n = 0.4, L = 10 
are shown in Fig. QJa) with thick lines. 

The behavior of the magnetic character of the impurity 
affects in turn all the magnetic properties in the system. 
Let us first examine the spin-spin correlations between 
the impurity site and the remaining sites along the con- 
duction chain. Results for the L = 11 chain with peri- 
odic boundary conditions (PBC), n = 0.364, are shown 
in Fig. EJa), (b) and (c) in the three regions where the 



spin of the ground state is S max , and S„ 



1, respec- 



tively. We also included in this figure results for L = 19, 
n = 0.421, and L = 25, n = 0.4, OBC, which show the 
same behavior. It should be noticed that in addition to 
the expected larger magnitude of these correlations in 
the S — region with respect to the ones in the other 
two regions, there are also qualitative changes. To detect 
these qualitative differences we have computed the static 
structure factor along the conduction chain: 



1,3 



(2) 



where j, I label the conduction chain sites, and q = 
(2tt/ L)n, n = 0, — 1 . Results where site j is 

restricted to the impurity site, that is just the Fourier 



FIG. 6: (Color online) Spin-spin correlations from the im- 
purity site along the conduction chain for several values of 
e in the regions where (a) S = S max , (b) 5 = 0, and (c) 
S = Smax — 1. The normalization (So So) = 1 was adopted, 
(d) Static structure factor along the conduction chain for var- 
ious values of e. Results for L = 11 chain, PBC, U — 4, 
n = 0.364, except otherwise stated. Results for L = 19 corre- 
spond to n = 0.421 and for L = 25 n = 0.4. 



transform of the correlations shown in Fig. [H^a), (b) 
and (c), are essentially the same as those obtained us- 
ing Eq. @. As it can be seen in Fig. Hfd), x(q) nas 
the typical shape of a FM order in the subspace of total 
S* — in the regions with S = S max and S max — 1, while 
it presents a peak at the smallest nonzero momentum in 
the region S = 0. As observed earlier, in Fig. [3jd), there 
are no traces of AF order in this region. 

Further information about changes in magnetic prop- 
erties caused by the impurity can be obtained by looking 
at the dynamical impurity susceptibility defined as: 

= £K*»I 5 W*°>I 2 *( W - ( E n - Eo)) (3) 

n 

where the notation is standard. Since the contribution 
from the remaining sites on the conduction chain is negli- 
gible, Si mp (ijj) is essentially equal to the total dynamical 
susceptibility once the contribution from localized spins 
has been subtracted. In all results below, the peaks have 
been broadened with a width 8 = 0.1. 

In Fig. EJa), (b) and (c) we show Si mp {uj) in the low-w 
part of the spectrum for various values of e in the re- 
gions S = S max , and S max - 1 respectively, for L = 11, 
n = 0.4 and L = 19, n = 0.421. Again, in addition 
to an expected difference in amplitude, a clear qualita- 
tively different behavior is noticeable. In this part of the 
spectrum Si mp (u>) presents small peaks at finite u> in the 
S = S max and S m ax — 1 regions while in the S = region, 
Si mp {u)) shows a strong peak close to u = 0. Results for 
L = 25, n = 0.4, e = —8 arc indistinguishable from those 
for L = 19. 

The distinct behavior caused by the impurity can also 
be detected by looking at the dynamical structure factor 



FIG. 7: (Color online) Dynamical impurity susceptibility 
Simp{u>) vs uj for U — 4, n — 0.4, and values of e in the 
region where (a) S = S m ax, (b) S = 0, and (c) S = S ma x - 1. 
Energies of the dominant peaks in the dynamical structure 
factor xil, 1 ^), (d) in the high- and (e) in the low-o; regions 
as defined in the text for various values of e. The size of the 
symbols is proportional to the intensity of each peak. Results 
for L = 11 chain except otherwise stated. 

on the conduction chain which is defined as: 

X (q,u>) = J2 K*nl^|*o>| 2 «(w - (E n - E )) (4) 

n 

where S* — (1/L) 5*| exp (ijq), and the sum extends 
over the conduction chain. Figure [7^d) and (e) show the 
energy and intensity of the centroid of the main peaks 
in the high- and Iow-oj parts of the spectrum for values 
of e in the three regions of total spin S above discussed. 
In the S — S max and S m ax — 1 regions the peaks with 
largest weight are those with energy close to Jh which 
corresponds to the magnon excitation. This behavior is 
strikingly different to the one in the S — region where 
the most weighted peak form a dispersionless band at 
the bottom of the spectrum. This band is reminiscent 
of the one found in gapped spin systems upon doping 
with nonmagnetic impurities. 14 In the present case, these 
low-energy peaks may correspond to magnetic excitations 
living in a "cloud" surrounding the impurity that can be 
observed in real-space in Fig. [DJb). 

V. DISCUSSION AND CONCLUSIONS 

We have applied both well-established and recently de- 
veloped numerical approaches to study the FKLM with 
an Anderson impurity in the FM phase on finite chains. 
We found that the magnetic or nonmagnetic character of 
the impurity is determined by a relationship between the 
impurity parameters and the Hund's rule exchange cou- 
pling of the manganite. As expected, transport occurs at 
the crossovers between the empty, half-filled or filled QD 
regions. 

The most important result of this article is the pres- 
ence of an intermediate singlet phase between the fully 
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saturated phases with S max (empty impurity) or S max — 1 
(doubly occupied impurity). This problem has some re- 
semblance with the problem of the existence of an inter- 
mediate phase between two ordered states in the frus- 
trated Heisenberg model on the square lattice. This in- 
termediate phase had been predicted by ED studies on 
small clusters^ The alternative view was a first order 
quantum phase transition between both ordered phases. 
Only recently was this controversy being settled favour- 
ing the existence of the intermediate stated but the na- 
ture of this phase is still a subject of active research^ 
Of course the physics involved in both problems is com- 
pletely different, but by analogy, in our problem, we 
could consider the possibility of a first-order transition 
between two FM states instead of the proposed inter- 
mediate nonmagnetic phase. By comparing the models 
involved in these two problems, FKLM and frustrated 
Heisenberg model, the former is much more difficult to 
analyze than the latter by numerical techniques since the 
size of the Hilbert space is much larger for a given cluster 
and moreover taking into account the convergence prob- 
lems discussed in Section [TTJ These difficulties prevent 
us to perform an extrapolation in order to decide if this 
intermediate phase really exists in the bulk limit or if it 
is just how the transition between the S ma x and S ma x — 1 
phases manifests in finite systems. 

In any case, the possibility of the existence of this in- 
termediate state is in principle interesting and important 
and it deserves further study. In this sense, there are 
three issues that should be considered. The first one is 
that this intermediate phase could be stabilized by some 
modifications of the model, for example by including the 
Heisenberg interaction between localized spins, or by re- 
placing the spin-1/2 localized spins by higher spin ones 
which, in addition, are also more realistic for manganites. 
The second issue is that even if this intermediate phase 
has a finite range around an impurity, a finite density 
of impurities could lead to a macroscopic feature. The 
situation here is analog to the presence of nonmagnetic 
impurities in the above mentioned gapped systems.— In 
these systems a single impurity attracts locally a spinon 
to the impurity and a finite density of impurities drives 
the system to a long-range AF order. These two issues 
are currently under study^ The third issue we would 
like to consider is related to the relevance of the present 



model to devices with finite dimensions. In these meso- 
scopic systems, as discussed in introductory textbooks, 24 
due to its finite size, many physical properties are differ- 
ent to that found in bulk systems. It is then relevant for 
these devices to capture short-range effects. 

Finally, we would like to provide a qualitative scenario 
to help understanding this nonmagnetic state. Let us as- 
sume that the system is in a low S z state. In the region 
where the impurity is empty or double occupied, both 
"leads" are relatively disconnected and each one would 
have a ferromagnetic state with spins polarized in one 
direction and the other with spins polarized in the oppo- 
site direction. Of course this state is degenerate with the 
one with reversed polarizations. Now, when the impu- 
rity is singly occupied, not only it would have a definite 
magnetic character but it would allow the crossing of one 
electron from one lead to the other where it would have 
then a "wrong" spin. This kind of magnons would then 
decrease the total energy both by increasing the kinetic 
energy and by decreasing the magnetic energy due to an 
effective AF interaction with the impurity. Of course, 
this gain in energy would not occur for the fully polar- 
ized system so it would drive the system to lower S and 
presumably to S m i n £^ It is interesting to notice that this 
scenario would then imply an enhancement of transport 
in the system which could be relevant for the devices 
mentioned in the Introduction. 

In summary, we present a prediction on the magnetic 
state of the FKLM doped with a magnetic impurity. This 
prediction could be experimentally verified on Cu-doped 
manganite nanotubes. These results could also be in 
principle reproduced experimentally on spin valves where 
manganites are used as ferromagnetic leads. We hope the 
present results will encourage theoretical studies to fur- 
ther characterize this proposed intermediate phase and 
to explore its presence in more realistic models for man- 
ganites. 
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